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INTRODUCTION AND SYNOPSIS 


The "nucleation in supersaturated systems" (such as the formation of 
fog in supersaturated water vapor, for example) was originally made amen- 
able to quantitative treatment by Volmer and Weber (ref. 1). To every 
saturation there corresponds a certain critical droplet size of the new 
phase of such a type, that the vapor is supersaturated only with respect 
to those droplets which are bigger than the critical droplet, but not to 
those which are smaller. The formation of fog is therefore contingent 
upon the origin of "kernels" or nuclei, i.e., droplets of precisely that 
critical size by a typical phenomenon of fluctuation. The frequency of 
such processes is, according to the relationship between entropy and prob- 

ability, proportional to e , where Acrit. "the energy required 

for the reversible creation of such droplet. Volmer 's treatment is 
briefly reviewed in section 1. The proportionality factor K, as yet 
indeterminate (in our equation (5))> was calculated by Farkas (ref. 2) for 
the case of droplet formation by a kinetic treatment, the results of which 
are fully confirmed (in section 2) by a more lucid method of calculation. 
The drawback of Farkas' calculations, as well as the arguments advanced 
by Stranski and Kaischew in connection with it (ref. 5)^ is that these 
writers ' first convert the elementary eqiiations of the kinetic theorem, 
each of which refers to the evaporation and condensation of a single mole- 
cule, in a differential equation which, when integrated, produce new and 
not always lucid constants. The change to the differential equation is 
risky because the ensuing functions of the molecule number n are at 
first def. led only for integral values of n and at the transition from 
n to n + 1 change frequently so much that the differential quotient 
loses its significance. By disregarding this risk Kaischew and Stranski 
obtained an incorrect result which differs from that of Farkas. On the 
other hand, the change into differential equation is entirely unnecessary 
(as will be shown in section 2). The algebraic equations for the indi- 
vidual processes give the wanted result by a simple, purely algebraic 
process of elimination. This method is shorter and less subject to errors 
than that of Farkas. Furthermore, the appearance of indeterminate 
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integration constants is completely avoided. Thus Farkas’ final for- 
mula, for example, still contains a constant which he himself designates 
as indeterminate, while in reality, an accurately estimable value cor- 
responds to it, which is in optimum agreement with the Volmer and Flood 
measurements . 

The next three sections deal with the origin of critical nuclei, to 
which the general thermodynamic analysis of section 1 is applicable as 
for the droplets.^ The first kinetic calculation of the thermodynami- 
cally indeterminate quantity K for crystals was made by Kaischew and 
Stranski (ref. This important investigation prompted the present 

study. With regard to the highly idealized crystal model, use is made 
of the simple cubic lattice, utilized by Kossel as well as by Stranski, 
which consists of nothing but cubic basic elements, which are in ener- 
getic reciprocal action only with its six neaxest neighbors. However, 
our results are largely independent of this special model conception. 

The kinetic analysis, like that of Stranski and Kaischew, results in a 
confirmation of Volmer 's formula. On top of that, we succeeded in 
defining the absolute value of K for this case too. 

Our algebraic method of eliminating the intermediate states not of 
direct interest affords an instructive representation suitable for 
the discussion of the particular nucleation process on the passage of 
an electric current through a network of wires of specific electric poten- 
tial differences at the ends of the network and given ohmic resistances 
of the individual wires forming the network. 2 

The whole discussion of the system of algebraic equations is then 
equivalent to an investigation of the conductivity properties of this 
network. This method produces in sections 4 and 5 comparatively simple 
and clear calculation of nucleation frequency for two- and three- 
dimensional nuclei. 


^Kossel's contrary opinion (Ann. d. Phys. (5), 21, p. 193^) 

stems from a misconstrued conception of the nature of thermodynamic 
considerations, which never refer to individual molecules but to those 
average values which in technically feasible experiments, come under 
observation. For example: the work of separation of the single molecules 

in a lattice plane may Jump back and forth arbitrarily; but in the evap- 
oration of the total lattice plane, only the mean separation work enters 
the balance of the thermodynamic process as heat of evaporation. 

^The possibility of such a representation was originally voiced by 
R. Landshoff in a conversation. Another, even more instructive repre- 
sentation is that of a diffusion process. (Cf. Volmer, Z. f. E., 55, 
p. 555, 1929.) But for the purposes of a quantitative treatment, our 
electrical pattern should be superior to the diffusion pattern, especially 
when a change from droplet to crystal is involved. 
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As an example for the application of the obtained results, the expla- 
nation and limits of validity of Ostwald's step rule are discussed in sec- 
tion 6^ Lastly (section 7)> mention is made of the unusual and rather 
general fact that in our electrical representation of the process of 
growth the resistances of all separate wires, which start from a specific 

+ 

state in the direction of growth, are given exact by Constant x e , 
where A is thermodynamic potential of this state with respect to the 
initial state (vapor, for instance). The kinetic interpretation of 
Volmer's formula (5) amounts then to indicating that the total resistance 
of the network is dependent solely on those pieces of wire which lie in 
the region of the point related to the critical droplet or crystal. 


1. THERMODYNAMICS OF NUCLEATION 


If n denotes the nmber of molecules contained in a droplet, F 
its surface and a its surface tension, the relationship between its 
vapor pressiire pn and that of a flat fluid surface (poo) reads 

dnkT In — = a dF (l) 

Poo 

where dn is the increase in the number of molecules corresponding to 
the surface increase dF. With the radius rjj of the droplet for spher- 
ical shape 

n = rn5 £ and F = 

5 m 

hence 


In Pn = 2^ i (2) 

Poo kTp 

is the critical droplet radius corresponding to the pressiore Pn- 

At given pressiare, droplets vith smaller radius evaporate, those with 
larger radius grow. A droplet which is exactly in equilibrium with a 
given pressure, according to equation (2), is hereinafter also designated 
as critical droplet or as nucleus corresponding to the particular pres- 
sure. A condensation of the supersaturated vapor can therefore take 
place only when a nucleus originates as a result of a fluctuation phe- 
nomenon associated with entropy decrease. 

According to the Boltzmann relationship between entropy and proba- 
bility, the probability for the appearance of such a droplet is 


4 


NACA OM 157^ 


_ 2 . 

proportional to e where S is the entropy decrease associated with 
the formation of a droplet of radius rn from a vapor of pressure Pn 
at constant volume and constant energy. If the manber of molecules con- 
tained in the vapor space is excessively great with respect to n, this 
entropy decrease is eqml to 1/T times the work A that must be per- 
formed in order to produce such a droplet in the vapor space isothermally 
and reversibly. This work can be determined, according to Volmer, by the 
following process : 

1. Removal of n molecules from the vapor space 

2. Ebqjansion of p^ to Poo 

5. Condensation on a flat fluid surface 

4. Formation of droplet from the fluid 

The sum of these four operations must give the wanted quantity A; 
but (l) and (2) compensate one another, which leaves 


A = -nkT In ^ + oF 
Pco 


Hence, with eqmtion (l) borne in mind 


A = ctf(i - ^ (5) 

\ F dn/ 

Since F = Constant x n^/^ it follows that S. ^ = 2. that is 

" F dn 3 


A = 





For the number of fog droplets produced per second, denoted hereafter 
by the letter J, we therefore expect 


- 

J = Ke 5kT 


(5) 


where Fn is the surface of the critical droplet corresponding to the 
given pressure p. The factor K still remains indeterminate in the 
thermodynamic study, and must be defined by kinetic analysis as origi- 
nally made by Farkas. The subsequently chosen method of computing K 
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is clearer from the methodical point of viev. Aside from that, the origin 
of crystal nuclei is to be treated also for which this equation ( 5 ) must, 
naturally, be applicable too. It will be seen that the factor K for 
fluid nuclei and crystal nuclei of equal order of magnitude is given by 
the gas kinetic collision factor. 

Regarding the differential quotient dF/dn in eqiaations (3) and (l), 
it should be noted that dF/dn is the mean growth of surface in the devel- 
opment of a molecule. For in the thermodynamic equation (l), dn 
still must always contain a multiplicity of molecules, althou^ the equa- 
tion is inapplicable as yet to single molecules. If this averaging of 
svirface growth per molecule is not carried out over a greater number of 
molecules, the surprising result is that the concept of vapor pressure 
loses its simple meaning for crystals, as shown by Kossel (ref. 5)^ 
because the increase of crystals in the growth of a molecule is mostly 
zero, but now and again very great too. 


2. FLUID NUCLEI 


Consider the following quasl-stationary condensation process. The 
vapor pressure p in a very large tank is kept constant by addition of 
single molecules . Droplets are then produced continuously which wovild 
increase infinitely without outside interference. To prevent this, each 
droplet, as soon as a certain number s of molecules is reached, shall 
be removed from the tank and counted. 

With regard to s it is simply stipulated that it shall be greater 
than the critical nvnnber n. The number of droplets per second counted 
under these conditions is termed "nucleation frequency. "5 In this pro- 
cedure, a steady distribution of droplets of various sizes will occixr 
within the tank, which must be examined a little closer. Suppose that 
Zy is the number of droplets containing exactly v moleciiles. The num- 
ber of free vapor molecules kept constant in our tank by addition is 
then Zi, while Zs is held to zero. If J is the mmiber of droplets 
counted per second, J may be regarded as a constant current that passes 
through all Z. 

Next, assume that: 

qydt is the probability that in time interval dt, one molecule will 
leave 1 cm2 of the surface of a drop of v molecules, aQdt on the 
- 

This term coined in the literature is somewhat misleading insofar 
as the actual number of nuclei formed per second is exactly twice as great 
because there is precisely a 50 percent probability for each nucleus to 
continue to grow or to evaporate. 
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other hand is the probability that one molecule from the vapor space 
condenses on a surface of 1 cm^, Fy is the surface of a droplet with 

V molecules, and Zy' = ZyFy is the total surface of all droplets with 

V molecoiles 

Applied to the constant current we get 


J = aoZy_x' - (lyZy* (for all v). 


Indicating 


^ (6) 

4.V 


the initial conditions read then 


Zy' = Zy_l' 3y - -SZ- Py (j) 

ao 

The factors 3 introduced by equation (6) increase monotonic with 
increasing v. For the critical molecule number v = n, Pn “ 
ry denotes the radius of the droplet with v molecules, then, by 
equation (2) 

p 2 om/L _ jJ\ 

rvj (g) 


The factor occurring in the exponent is indicated by 


= 2oM 
pRT 


(8a) 


^Itaf ortunately, the notation Zy ' and Zy in the Stranski and 
Kaischew article are enterchanged relative to Farkas' report. ¥e follow 
Farkas' notation. 
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In order to eliminate from the equations 


Zv+l' = Zv'Pv+1 - ^ 3v+l 

^v+2' = 2 >v+i'Pv+2 - ^ 3v+2 

> 


Zs' = Zs-l'Ps - 4 


(7a) 


the factors Zy+lS ^v+2' f • • • Zg.i’, the first is divided by pv+1, the 
second by Py+lPv+2^ etc., the last one by Py+lPv+2 • • • 3s* When all 
the thus obtained equations are added up, all the Z' values lying between 
Zy ' and Zg ' cancel out leaving 


Zs* 

3v+i3v+ 2 • • • 3s 



3v+i3v+2 


+ 


Pv+i3v+2 


3s-l, 


With it the nucleation frequency J is known as soon as one of the values 
of Z' is given. In view of the calculations for the crystal nucleus this 
method of solution is somewhat modified as follows; Through the multi- 
plications equation (7a) takes the form 


'^i+1 = - JRi 


(9) 


with 


<1.1 


Zj' 

P235 . . . 3 i 


and 


Ri 


1 

^0^233 • • • 3 i 


(9a) 


The quantity 3.^ arises from the corresponding Z' values by divi- 
sion by the product of all the 3 values which occur during the succes- 
sive growth of the droplet characterized by subscript i from single 
molecules. (By this method the equations are divided by the common fac- 
tor ^ 2^3 * * • 3y*) style of writing (equation (9)) of the equation 
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system indicates that the current J flows from point i toward point 
i + 1 under the influence of the voltage difference by over- 
coming the ohmic resistance Visualizing a series connection of 

resistances R]_, B-2f etc., the entire nucleation current J can be 
regarded as a current driven by a given potential difference throu^ this 
chain 


J(Rv + Rv+1 + • . • + Rs+l) = ~ ^s 


Now 0]_ is directly equal to Zj_' and Og equal to zero. The whole 
problem therefore consists in adding the separate partial resistances. 
Now it is seen that the individual Pv values increase in such a way 
that exactly equal to unity, while the preceding ones are all 

smaller and those that follow all greater than unity. Up to the value 
Rji the partial resistances consist, therefore, of a product of integral 
factors which are greater than unity; on above R^ the additive factors 
appearing are all less than \xnity. As a result the R^ values plotted 
against i have a distinct maximum at i = n. Owing to the importance 
(8) of the quantities p the exact term for a partial resistance R^ 
reads 



The sum of the reciprocal radii occurring here in the exponent is 
replaced by an integral with respect to the quantity 


Xv 



( 10 ) 


The integration variable x indicates^ therefore, the ratio of a 
particular droplet radius to the critical radius . By solution of equa- 
tion (lO) vith respect to v 

V = n(xv)^, dv = 3nx^dx 

hence 




dv _ 3n 


v=l 


■n 


r ^ xdx = I ^(xi^ - xi^) 

Jxp ^ 
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In addition 


i - 1 = n(xi^ _ x^P) 


Indicating for abbreviation 


A' 


on 



5kT 


the term for partial resistance Rj^ reads 

Ri = i e^' {(3xi2-2xi5)-(5xi2_2x;L^)J' 
ao 


( 11 ) 


Replacing the summation over the partial resistances also by an 
integration, leaves 

r Rydv = 5n f R(x)x^dx = ^ e“^' (3xi^-2x]_^) F gA' (5x2_2x^)x2(jx 
J 1 J ao J xi 

At X = 1 the integrand has a steep maximim of the order of e^' . There- 
fore we put x = l+|, i.e., 5x^ - 2x^ = 1 - 51^ - 2|^, and get the 
integral 

J e-A’ (5|2+2|5) (1 


The variation of the integrand is represented in figure 1. The 
factor A' is fairly hi^, say about equal to 20 to 50 , in practical 
cases, as will be shown later. So, without appreciable error the above 
integral can be replaced by 


e-5A'|2 

Then, the total resistance (5x]_2 _ 2x3_5 compared to unity being disre- 
garded in the exponent) reads 




( 12 ) 
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With this the thermodynamically obtained expression for the nucleation 
frequency of the indeterminate constant K is defined. 

The final result is 


J = 




(15) 


Against this calculation the objection mi^t be raised that the 
formula (l) had been applied to droplets of as low as two or three mole- 
cules, for which the concept of surface tension is certainly perfectly 
meaningless. But, when considering the curve of the partial resistances 
in figure!, it is clear that the resultant total resistance is definitely 
defined by the partial resistance in the neighborhood of v = n. There- 
fore it is practically immaterial whether the partial resistances at the 
start of the chain had been chosen by a factor 100 too great or too small. 
Equation (13) is exactly identical with Farkas’ formula (ref. 2), when 
bearing in mind that his constant C on the basis of its introduction 
(p. 239) has the significance Since Farkas did not notice that 

the extrapolation of his formula to droplets of only two or three mole- 
cules is positively unobjectionable, he failed to recognize the signifi- 
cance of this constant. 


In comparison, the calculation of Kaischew and Stranski (ref. 3) does 

. , dZv ' 

not seem to be entirely acceptable. They replace Zy_]_ - Zy by 


dv 


which serves no useful piirpose in the subsequent calculation, since no 
integration along this differential quotient is ever made. It merely 
obscures the significance of their constant C which simply is -Zp'. 


But, contrary to Farkas, they use the calculating method of logarithms 
and subsequent substitution of the differential quotient for the differ- 
ence quotient for great v also. This certainly is inadmissible in 

the range of small , where the logarithmic term changes rather con- 

dZv ' 

siderable even at minor changes in . 

dv 


The formula obtained for J is now compared with the Volmer -Flood 
measurements on fog formation at adiabatic expansion of water vapor. 

The factor Zp' is, by assumption, eqml to the total surface of the 
free molecules; ®’0^l' signifies thus twice the mmiber of gas kinetic 
collisions per second between the Zp vapor molecules . From the mean 
free path length 1. and the mean molecular velocity v the number of 
collisions per cm^ of vapor space follows at 
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Since is inversely proportional to the concentration N, we get 


aoZi' 


n 2 -JL- 
IoNq 


L\ I— 

rtR PqIq \flA 1\fl 


= 5 X 10' 


,22 _1 Pf 

ZqVM 1 '^ 


where 

N number of vapor molecules per cm5 

L Loschmidt number 

p vapor press\jre mm Hg 

IjQ free path length at 0° C and standard pressure 

For the nmber of molecules in the critical droplet we get 


whereby In 



X. 


n = — r„^ — 


M 


L 


4it /2aM l\3 p ^ 

3 \pRT xJ M 


For water (p = 1 g/cm3; o = 75 dyn/cm) 



For the attainable super saturations (x «= I. 5 ) n amoxants to about 
100 molecules . 
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Posting 


A ' - i fis. = L2/'2^ l-Y 

5 kT 5 RTVpRT xj 

for A' the formula for J reads 

2 

In J = 1*9 + In ^ + 2x t 2 In X - (13a) 

IdM5/2<,5/2 ==2 

Poo in mm Hg 

This result is then ccanpared with the Volmer -Flood measurements on 
water at temperatures T of 260° and 275°. All measurements at 
T = 270°, Poo = ^ mm Hg, Z_ = 10“5 cm, cr = 75 dyn/cm are entered below 
the logarithm. Hence 


In J = 52.5 + 2x + 2 In X - 5-7^ X 

The curves obtained for In J are shown plotted against x in 
figure 2 for T]_ = 275-2° and T2 = 26l.0°. But there is a certain 

uncertainty as to which value of J is to be designated as condensation. 
According to the graph the curves intersect the x axis at such a slope 
that it is practically immaterial, when defining the critical supersatu- 
ration, whether J = l(ln J = 0) or J = 10 (in J = 2.3) is plain fog. 
Choice of the intersection- point of the curves with the straight line 
In J = 1, gives the following values for the critical supersaturation, 
which can be compared with the measurements 






Pn 



T 

c(dyn/cm) 

X 

calculated 

Poo 





1 

calculated 

measured 

Curve 1 

275.2 

75.25 

1.46 

4.30 

4.21 

Curve 2 

261.0 

77.28 

1.64 

5.14 

5.03 


A' amounts to 55 to 56 


Worthy of note is the insensitivity of the theoretically computed factors 
to errors in the calculation of b-qZ^^.' • Even a factor 10 would change 

the constant 52.4 only by 2, i.e., practically no change at all in the 
result. 
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Since all further measurements on other substances in the Volmer- 
Flood report have been compared with those measurements on water and 
gave a good confirmation of Parkas' formiHa, it is concluded that the 
present formula (13) reproduces the whole available test material very 
satisfactorily . 


3 . THE LINEAR CHAIN 


Preparatory to the problems of the actual crystal growth, the fol- 
lowing process is analyzed: It is assumed that the rectangular area ABCD 

is the base of a simple cubic crystal, on which as the beginning of a 
new sxirface, a layer of edge lengths z and 2 is available and on which 
the (2 + l)th chain of length z is included in the growth. The growth 
of this new chain is analyzed. Figure 3 represents the stage in which 
exactly k = 3 atoms of the (2 + l)th chain are condensed. The diffi- 
culty of forming nuclei here is due to the fact that during the start of 
a new chain the first and possibly also the second and third atom are 
less solidly bound than those following, which are all boimd with the 
same energy (repeatable steps, according to Kossel, bond at "half crystal" 
according to Stranski). So, unless there is too much supersaturation 
after a ccanplete chain has formed, there is a considerable lapse of time 
before - as the start of a new chain - a linear nucleus capable of growing 
has formed. The energies, with which the single atoms are bound in the 
successive formation of the chain, are indicated with <pi, cp2> • • • 

9k> .... Then the possibility q}tdt that the k-th atom evaporates 
as a result of the thermal motion in time interval dt on a chain con- 
sisting of k atoms, is given by 


. ?k 

qij = F(T)e‘ kT 


(14) 


On the other hand, the possibility adt, that a further atom settles on 
the chain, is independent of k and solely given by the external vapor 
pressure. It is assumed that there is no slip of atoms at the crystal 
surface. In that event a is essentially equal to the number of vapor 
atoms per second arriving at the surface of a single crystal atom. The 
quantity a introduced here follows from the uq (of section 2) by 
multiplication with the surface atom. 

We put 


_ i_ 

a = F(T)e 


(15) 
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hence, where the energy \tr is the measure of the external vapor 
pressure. 5 The whole mechanism of growth is governed by the factors 

i . e M ( 16 ) 

which for normal growth at the beginning of the chain (k = 1) are sub- 
stantially smaller and for greater k a little above unity- In conjunc- 
tion with Stranski and Kossel, this behavior is then schematized so that 
Pi is regarded as very small compared to unity and all other p values 
as equivalent and greater than unity- For the investigation of the growth 
of a chain the following steady process is analyzed: A space under con- 
stant vapor pressure contains a very large number of crystals which are 
in the stage of growth represented in figure 3- But the new chain in 
the process of formation may be of any possible length and assume any 
possible position on the raised side of the rectangle. The number of 
crystals on which the new chain has exactly the length k and is at a 
specific position at the growing edge is indicated with n^; correspond- 
ingly, the number of crystals arising from the crystals of the type nj^ 
due to deposition of an atom at a certain end of the chain k, is indi- 
cated by njj-^p. 

By partial current J’ is meant the excess of the growth process 
per second which lead from the nj^ crystals to those of the type 
through the evaporation processes, which lead from to n^- For 

this specific partial current 

J’ = njjH - 


Each chain has then two possibilities of adding an atom corresponding 
to its two free ends. In the two positions of chain k in which one end 

coincides with one end of the base, there is only one possibility of 

build-up. Since, for the chain k, there are (z - k + 1) various posi- 
tions possible, there are altogether 2(z - k + 1) - 2 = 2(z - k) partial 

currents J', which collectively lead from all the crystals with chains 
of length k to those with chains of length k + 1. But in the case 
of transition from nQ to np there are only z partial cixrrents 
corresponding to the z deposition possibilities of the first atom of 
the new chain. This branching of the current is represented in figure 4 
for z = 6. 

Now the not entirely exact assumption is made that all partial cur- 
rents leading from k toward k + 1 are eq^ui valent. Since their sum 

^At absolute zero point the heat of vaporization would have to be 
used. As a rule, signifies a thermodynamic potential. 
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gives the total current J, each is eqml to 


J- = — J 

2(z - k) 


It is readily apparent from figure 4 that this assumption cannot he 
rigorously correct on account of the equations of continuity between the 
partial currents. Owing to equation (I 7 ) this assumption corresponds to 
the assmption that all positions of the chain k are identically 
frequent . 


Thus, on this premise the steady state is described by the equations 




ni = noPi 

_ J £1 

~ a z 

n2 = n^pg 

J 32 

2a z - 1 


J ^k+1 

“k+1 - "k3k+l - 7TT 


( 18 ) 


~ ^z-l^z ” 2a 

These equations are treated the same way as those of the droplet forma- 
tion (in section 2 ), by regarding them as equations for the passage of 
current through a series of specified partial resistances, ^y division 
of the k-th equation by the produce Pk+1 = 3l32 • • • 3k+l# they take 
the form 


«’k+l = ^k - J8k 


(19) 


where the individual potentials and partial resistances indicate 
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d-v = 




^1^2 • • • 3 


Rk = 


2a(z - k)P2^ . . . 

but *0 = 




'*0 = i 


(19a) 


Thus, in figiire 4 the electric potential of a specific state is 
represented by the quotient of the nmber nk of crystals in state k 
and the product Pk = PiP2 • • * Pk P values of all atoms bound 

in this state. Specific experimental interogatory forms are synonymous 
with the corresponding statements regarding the electric potential dif- 
ference placed at the ends. However, it is to be noted that, in contrast 
to the electrical picture, the absolute value of the potential in the 
growth process itself has a well defined meaning 

•tk ^ 

— and respectively 

I^k ^k-1 

namely are the number of Individual processes taking place in unit time 
from k to k + 1 and k to k - 1. 

As application of (l9) the actual linear nucleation as well as the 
growth of the rectangle layer about a whole chain is now analyzed. 


A. Linear Nucleation 


The procedure for defining the nucleation frequency is the same as 
for the droplet formation. All the crystals for which the chain has 
reached a certain arbitrarily chosen length s are removed; s is to 
be very small compared to the length z of the edge. The number of 
crystals removed per second is called linear nucleation frequency. 

Hence we put = 0 and find 


J = 


R 


with Oq - riQ 


and 

R = + . . . + ^s-1 

J^(k + 1 + 1 + . , _ + 1 

2a\z (z - 1)3]^ (z - 2)3 iP2 (z - s + l)pi . . . 3s-l 
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With the specialization P2 = Pj = . . • = P and. because s « z 

R = — 12 + — + ps-3 + . . . + 1) 

2az L J 

(3 - l)3j^3“'V 

is applicable also. 

Disregarding the 1 next to p^"l and the 2 next to l/Pi^ leaves 
1 P 

R = — , hence the frequency of the linear nucleation at one 

2 az (p - l)pi 

of the Hq edges, independent of s 

— = 2 az Pi ( 20 ) 

^0 3 

The factor 2 ^ ~ ^ Pi, small compared to unity, is regarded as a 

probability that one of the atoms striking the edge (their number per 
second amounting to az) grows up to a new chain. 


B. Deposition of a Whole Chain of Length z 

In this event all the partial resistances from Rq to Rz_i must 
be added up 

^ ^1 ^z-1^ “ “ ^z 

or 

\ _ >^z 

1 = no 

• 3z-l/ Pi32 • • • 3z 

The first partial resistance 2/z in comparison with — 

(z - l)Pi 

can always be disregarded. Putting ^2 ~ • • • = Pz ~ Set with 

the abbreviation 


2 + 


2 aU (z - l)3i (z - 2 )Pip 2 


P1P2 
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nz = UqPiP' 


J s 

2a 


The sum S 2 does not lend itself to elementsiry evaluation. The 
approximate value 


Sz(p) = 


z In p 


(21a) 


used in the following is obtained by the following consideration: 
Replacing the s\jm (21) by an integral gives 


hence, with the substitution x = z - 


z In p 


In p(z-l) 


y 

z In p 


The approximate value is obtained by disregarding which is 

z In p 

small compared to unity, which, however, presents only a rough approxi- 
mation near the upper limit of the integral. Equation (22) enables the 
deposition of a whole chain z to be treated as an elementary process. 
The equation ( 17 a), valid for the actual elementary process, is simply 
replaced by the relation 


J — nQAz ” 
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whereby 


A 2 - 


2apip 


z-1 


= 2az In p 


31 


Bz = -S 2az In pP"Z 


(25a)6 


Az is slightly dependent on z, while Bz decreases exponentially with 
z. Both quantities become equivalent at a critical value of z, which 
is denoted by m, and is defined by 


PlP“-l = 1 or ^ = p“i (24) 
3l 

m is that chain length which is precisely in equilibrium with the exter- 
nal pressure. According to ( 16 ) the definition (24) of m is equivalent 
to 


91 + (m - 1)9 = mitr 

or (24a) 

9 - 'If = - (9 - 9l) 
m 

The mean evaporation energy of the "critical chain m" is equal to 
the energy t characterizing the external vapor pressure. 


4. TWO-DIMENSIONAL NIICLEUS 


Equation (22) makes it possible to analyze a chain of length z as 
an element, through whose deposition or evaporation the growth of plane 
nuclei or of whole rectangular plates is controlled. In this instance 
the growth of a plane nucleus on a given base of edge lengths i and 
k is involved. A specific stage of this growth is represented in fig- 
ure 5 . The bonding energy of a single atom on the smooth base ("bond to 
one neighbor") is denoted with 9 qj 9i and 9 have the same meaning 
as in section 5. Accordingly, there are 

To-'tf ^1-^ (D-i 

PO = e 4T ^ Pi = e kT , P = e kT (25) 

The energy required to detach the whole plate (i, k) from the base 
is then 

90 + (i + k - 2)91 + (i - l)(k - 1)9 


°In this calculation it is assumed that at no time two nuclei are 
simultaneously existent on the same chain and then grow together to one 
chain. When z is not extremely great, this assumption is well justified. 
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From the assumption that P 2 = Pj = • • • = 3 are all equivalent 

and independent of the position of the deposited molecule on the base, 
inevitably follows the condition 


cpo + 9 = 2cpi 


hence also 


3q3 = 3i2 


( 25 a) 


The total bonding energy of a structure must be independent of the 
manner in which the growth takes place. Applied specifically to a system 
of 5 atoms on the base as in figure 5^ the binding energy for growth in 
the order of 1, 2, 3, is cpQ + cpQ + qp, but for the sequence 1, 2, it 

is CpQ + cp 2 _ + cpQ_. The equality yields the above relationship. The evap- 
oration energy of the whole plate (i, k) is therefore 

ikp - (i + k)(cp - cpi) 


Visulize a column of cross section i x k consisting of I whole 
atom layers, the deposit of the {l + l)th layer being located on a 
rectangle s x z. In analyzing the f\ill growth of this deposit into a 
whole layer the procedure is the same as in section 3- A multiplicity 
of columns and rectangles of every possible size and position is assumed, 
with ng ^2 denoting the number of those at which the nucleus (s,z) has 
a specific position on the base. The number of the possible positions 
is, obviously (i - s + 1) (k - z + 1) . The total number Zg ^2 4he 
columns with a plate (s,z) woiild then be ng^ 2 (i - s + 1) (k - z + 1), if 
it is assumed, as in section 3, that each position of the rectangle (s,z) 
on the base occurs with the same frequency. For the current Js^z^ 
leading from s,z to s + l,z, there are altogether 2(i - s)(k - z + 1) 
possibilities, namely, two each for each specific position of the rec- 
tangle (s,z), with exception of those positions at which it lies to 
the left or right at the edge. In these cases there is only one possi- 
bility for depositing a new chain. It is assumed again that these 
2(i - s)(k - z + 1) partial currents, all of which lead from s,z to 
s + l,z, are equivalent. 


Now, in order to describe the growth of the plane crystal after the 
foregoing argments with the above eqmtion (22), n^ 

no by Ug^2 J ^ 


^s+1, z> 


2(i - s) (k - z + 1) ’ 


is replaced by 
while Jar,' is 
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to denote the ctirrent from (s, z) to (s, z + 1 ) . In the steady state the 
plane growth is then governed by the equations 


ns+1,2 = - 4^(1 _ 3)(k . ^ + 1) 


s. z 


s,z = 1, 2, . . 


^s,z+l ” ^s,zPiP 


On the other hand 


s-1 _ 


Js, z ' 


4 a(i - s + l)(k - z) 


(26) 




is applicable. 


nil = nooPo - ^ PO 


(26a) 


Now the content of these equations is described by a discussion in 
the s-z plane (fig. 6): 


Suppose that a certain lattice point s,z represents the crystallite 
defined by the edges s and z. Thus, in figure 6, for example, the 
point A corresponds to the crystallites 3 x 2 . The current J„ _ 

O ^ u 

flows then horizontally from s,z toward s + 1, z, Jg 2', but vertically 
upward from s,z toward s,z + 1. The whole lattice extends to s = i 
and z = k. The problem then consists in computing a total current J 
that enters at (O.O) and branches off in partial currents Js,z 
Js,z’> according to equation (26). The obvious method is to regard the 
entire figure 6 as the image of a material network through which passes 
a current J under the effect of a certain electrical direct current 
voltage. To complete the picture, the several equations of (26) must be 
expressed in the form of Ohm's law 


^s+1, z - < 3 >s,z " Js,z^s,z 
®s,z+l = 5 >s,z “ Js,z*^s,z* 


(27) 


which describes the current in the separate pieces of wire of the network 
in figure 6. It can be accomplished by dividing the first equation of 
(26) by the product P of all the p values occvirring in the build-up 
of the plane (s + l,z), that is. 


Ps+l,z = 3oPi(s+z-l)ps(z-l) 
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The potential <I>s,z the partial resistances Rs^z and Rs,z' 

become then 




n 


s , z 


p^p^(s+z-2)p(s-l) (z-1) 


■} ^ 0,0 = no,0 (27a) 


R 


4a(i - s)(k - z + 1) Pg+i^z^ 




Rs,z' = 


' 4a(i - s + l)(k - z) Pg^z+l 


(2Tb) 


The whole problem now consists in the discussion of the electrical 
properties of the network built up of partial resistances (27b). 

Introduction of the approximate value (21a) for S2 gives 


*^S Z 

’ 4a(i - s)(k - z + l)z In p p^p^s+z-lpsz-s-z 

Introduction of the critical edge length m defined by equation (24) 
and with the relation ^ ^ following from equation (25a) , the 

second factor of the above Rs,z becomes 


pmpm(s+z) -sz 


_ pm+m2 -(s-m) (z-m) 


With a system of axes turned through 45° (cr along the diagonal, C at 
right angle to it) 


s 


z 


a + U 


( 28 ) 


becomes 


(s - m)(z - m) = (a - m)^ - 
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hence the factor of Rs,z question 

pm+m2p-(a-m)2p^2 


The relationship of ^ implies: if the partial resistances R„ „ 

which are met by a normal to the diagonal, are examined, they are seen 
to have a sharp minimum on the diagonal itself (at C = O) . One, two, 
three lattice points away from the diagonal, the resistances increase by 

the factors p, p^, p^, etc. Thus practically only the diagonal s ^ z 

of this lattice is conductive, while along the diagonal (change of a 
at ^ = 0) , the resistance has the same sharp maximum at a = m, i.e., 
at the critical edge length. There the partial resistances drop with 
the distance from this point to the p-1, p-^, p-9-th fraction, when 

the distance from a = m amovints to one, two, three lattice points. 
Therefore, the entire current must flow along the diagonal in a narrow 
gorge, characterized by the factor which in turn leads over a very 

steep "pass" at a = m. 

The over -all resistance R of the entire network is practically 
concentrated on the resistance 


R Rjn m = (29) 

^ 4a(i - m) (k - m + l)m In p 

at the height of the pass. Naturally, it is assumed then that the spot 
s = z = ra still lies substantially within the lattice. A further discus- 
sion of the branching conditions in the immediate vicinity of the spot 
m,m could contribute to the expression R no more than a factor close 
to unity, which, however, would be useless for the purposes involved here. 
Now the two questions concerning the frequency of plane nucleations, as 
well as the growth of the rectangular col\min i x k around a whole layer, 
can be answered exactly as for the chain in section 5. 


A. The Formation of Plane Nuclei on a Rectangle Base 

To determine the nucleation frequency, all the columns for which 
the sum of the edges s + z has reached an arbitrarily specified value 
s + z = n are withdrawn frcan the vapor space and counted, that is, the 
lattice points lying on the straight line s + z = n are grounded. The 
saddle point s = z = m is found to be located still within the thus 
out-off triangle and also that n and hence m itself should be very 
small with respect to the edges i and k of the base. The result is 

the satisfactory approximation: R = ; hence, the nucleation 

Aaikm In P 
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cmrent for a single one of the nQ crystals is 

— = Aailm In (30) 

"o 

The two factors 3 ™ still remaining axe tied to one another 

q>-q>i 

by the relation 3°^ = ® t according to eqimtion (24), with the aid of 
which one of the two can then be eliminated from equation (30) . 

Elimination of 3 leaves 

m(T-9l) 

^ = 4aik ^ e kT g kT (30a) 

^0 kT 

Here, as it should be. 


half the free -edge energy 
kT 


is in the first e power, because ^(<p - cp^.) free-edge energy 

per atom, 4mi((p - cpi) is, therefore, the total free-edge energy. The 
2 

plane energy accompanying the formation of the chain (at both of its 
ends!) is in the second e power. The factor before the already thermo- 
dynamically required e function is extremely simple: its order of 

magnitude is defined by the number aik of the vapor atoms per second 

4(cp - cpi) 

impinging upon the plate ik. The then still remaining factor 

kT 

is nearly equal to j-op _ h^at ^ significance for the 

RT 

only interesting order of magnitude of J. 

On the other hand, when m substitutes for m + 1, the elimination 
of m from equation (50) gives 


(<P-Tl)^ 

J = 4aik — e (kT)^ In 3 (50b) 

kT 

This eqixation gives the dependence of current J on the experimen- 
tally directly defined supersaturation 3- 
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B. At the formation of a whole plane i X k it is necessary to ana- 
lyze the entire rectangular network on which the current J is introduced 
at point (0,0) and channeled off at (i,k) 


*^0,0 - “^iik = JRi,k 


where 


is the over -all resistance of the entire rectangular network. 


Here, also, the only ca£e of interest is that where the point (m,m) 
lies still within the net; Ri^k> therefore, according to equation (29)^ 
may be replaced by Rm,m- If use is made of the relations PqP = Pi 2 


and -^ = 

Pi 


in the product Ps,z> 


s , z 


then becomes 


Ps,z = PoPi®'^^"2p(s-l)(z-l) = p(s-m)(z-m)-m2 


Thus, by equations (27) and (29) 


Hi,. = - J -7- g^r’V 


4 a(i - m) (k - m)m In p 


( 31 ) 


On the basis of this equation ( 31 )^ the accumulation of a whole 
plane can henceforth be treated as elementary process. The equation 
applicable to it reads 


^ = ^0,0^i,k - ’^i,k®i,k 


(32) 


where 


Ai k = W(i - m) (k - m)m In pp-m-m^ 


(32a)''' 


Both qiiantities are equivalent at - (i - m) (k - m) =0 or 
ik = (i + k)m. For square plates, which are practically the only ones 


7 

Multiple nucleation is excluded again. At very great i and k 
the foregoing is therefore no longer applicable. 
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occiirring at the subsequent growth of the spatial crystal, it results in 
i = k = 2m. The probability for evaporation and condensation will not 
be identically great until the edge length of the plane is twice as great 
as the critical chain length m. Naturally, this result was to be fore- 
seen by reason of the fact that the equilibrium vapor pressure is defined 
by the mean evaporation energy, as already predicted by Stranski and 
Kaischew. It is to be noted that Ai^k is very slightly dependent on 
i and k, while decreases rapidly with increasing size of the 

plane. 


5. THE CRYSTAL NUCLEUS 


After these preparations the quantitative treatment of the nucleation 
frequency for spatial crystals is easy. Again visualize in a vapor space 
a large number of box-like crystals of all possible edge lengths i, k, 

I in steady distribution so that the vapor pressure remains constant 
and that all crystals, as soon as they have reached a certain size, are 
removed from the space and counted. The number of crystals with the edge 
lengths i, k, I is assumed at k Z* They may, for example, change 

to crystals (i + 1, k, Z) by gathering of a plane (k, Z), that is, this 
plane can be deposited on two different sides of the little box. The law 
for this process was defined in equation (51). Replacing nj^ k 

Zi+l,k,Z; nQ,0 ^i,k,Z> ^ results 

in 


^i+l,k,Z - ^i,k,ZP 


(k-m) (Z-m) -m^ 


J-i V 7 3 “+ (k-m) (Z-m) 

(55) 

8 a(k - m)(Z - m)m In 3 


where Ji^k,Z indicates the partial current that leads from i>k, Z to 

i + l,k, Z. By the method previously used several times, equation (55) 
can again be put in the form 


‘^i+l,k,Z 


'^i,k, Z “ '^i,k, Z^i,k, Z 


Oh) 


and identify it as the electrical current in a spatial network whose lat- 
tice points have the potential ‘I’i^k, Z ^ which Bi k Z ohmic 

resistance of the piece of wire that leads from i,k, Z to i + l,k, Z. 

The transition from (55) to Oh) is accomplished again by division of (55) 
by the product Bf+i^k, Z ^T1 the 3 values which occur diiring the suc- 
cessive development of the state (i + l,k,Z) from (i + l)k,Z single atoms. 
The factor for Zi,k,Z equation ( 53 ) is precisely the product Pk,z 
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of the 3 values for that plane (k, Z) which is newly added in the par- 
ticular reaction. The product \ can be constructed as follows: 

On a single atom free columns of i - 1, k - 1, Z - 1 atoms are depos- 
ited along the coordinate axes, each of which gives the factor pQ. The 

spread-out rectangle sides are then filled out. It yields (i - 1) (k - 1) + 
(k - l)(Z - 1) + (Z - l)(i - 1 ) times the factor and, in addition, 

(i - 1 ) (k - l)(Z - 1 ) times the factor p. Hence, altogether 


Pi,k,Z = Poi+k+Z-5p^(i-l) (k-l)+(k-l) (Z-1)+(Z-1) (i-l)p(i-l) (k-1) (Z-1) 


Po Po m 

for which the relations — = — = p“ give 

Pi P 


and 


p , _ p-m(ik+kZ+Zi)+ikZ+3m-l 
l,k, 6 


By division with this quantity, we obtain in equation (34) as factor 
of the final term for the partial resistance 


Ri,k,z = 


8a(k - m)(Z - m)m In p 


pm(ik+kZ+Zi) -ikZ+m'^-2m+! 


(55) 


To clarify the behavior of the exponent visxoalize a perpendicular 
line dropped from the point i^h, Z of the network on the space diagonal, 
OyOyO to indicate the foot of this perpendicular. 

Putting 


= a -h rn 


= a + r^ 


= a + r 


3 


(56) 


the construction given then 


r;L + T2 + r3 = 0 


( 36 a) 


28 


NACA OM 1574 


along with 


+ ’^2^3 + Vl 


- i (r 2 + r,2 t r,2) 




^ (36) and (36a) 


ifc + kl + Zi = 3 (j 2 - |r2 


and 


ikl = 0^ _ ~or2 + 


With this the exponent of equation (55) reads 


m(ik + kZ + Zi) - ikZ = 3cr^ni - 0 ^ + ^(0 - m)r^ - 


■1^2^3 


Disregarding the practically nonessential term rj_r 2 r^ (the sur- 
rounding of the diagonals being considered), the conditions for Ri^k, Z 
are the same as before in section 4 for the plane lattice. The factor 


^"(0 m^r^ 

p2'‘ ~ ' in the region 0 > m solely considered here, effects such a 

rapid rise of the resistance on leaving the diagonal, that the current 
can flow practically only on the diagonal r = 0. On the diagonal itself 
the factors p3cr2m-03 has such an enormously steep maximum at a = 2m, 
that the entire voltage drop along the diagonal is practically defined 
solely by the partial resistance 


I^2m,2m - 


8am^ In p 


54m^+m^-2m+l 


Owing to 'I>i 1 1 = 1 1 nucleation frequency is therefore 

defined at 


J = 8aZi^i^im5 In pp-4m3-m2+2m-l 


(37) 


Elimination of p by means of 


pm = ^ = e > m In p = ^ 

Pi ' ^ kT 
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leaves 


J = 


2 9 - 9i 

kT 


l^m2 (cp-(pj_) 


kT 


kT 


9 ^ 1 / 1 

" m 


(37a) 


The factors deciding the order of magnitude of J are aZ]_^2.,l 
the first of the three e-functions; aZ]_^]_^3_ is essentially (like the 

factor aQZ-j^' in equation (15) for droplet formation) the number of gas 
kinetic collisions per second. The first e-power is synonymous with the 

_ 

factor e 5kT of the thermodynamic formula (5). In fact, — (cp - cpi) 
is the surface energy per atom; the total surface energy of the cube of 
critical edge length 2m, therefore, is equal to (2m) 2 x 6 x ^(q> - q>i) = 

12m^((p - 9j_) > "the third portion of it stands, as it should be, in the 

exponent. The exponent of the second e-function indicates, as shown in 

1 Gd.ff0 6ii6r^^ 

section 4, = — for a critical plane nucleus. This factor 

2 ' kT 

occurs in similar manner in the report by Stranski and Kaischew too. 
Admittedly, its appearance hinges on the exact knowledge of the factor 
m, as is apparent from the fact that in equation (37) the term with m2 

can be made to disappear completely, if m is replaced by m - For 

the problem involving the critical supersaturation the second and third 
e-functions are ignored. 


6. THE OSTWALD LAW OF STAGES 


This law states that in the formation of nuclei from supersaturated 
vapor the liquid phase is separated first, as a rule, even when the tem- 
perature of the vapor is considerably below the freezing point. Our 
results on the nucleation of liquid (13) and solid (37) nuclei enable a 
theoretical foundation and a quantitative improvement of this law to be 
made. 


Omitting the last two e-functions in (37a) and introducing the rela 
tion for Volmer*s exponent of equation (5) 

A" = ^ = ^°^^(9 - 9i) 

3kT kT 



30 


NACA OM 157^+ 


we get by equation (3Ta) on the crystal nucleus 


Crystal: Jcrystal = 2aZi^l^lA"e-A 


while equation (13) produced 


Droplet. : J, 


droplet 


ap^l'' 


.-A' 


n 


The factor Z-j_' was, according to section 2, the number of vapor 

molecules multiplied by its surface. Since, according to section 3> a 
arises from Eq by multiplication with the atom surface, aZ]_^l^l and 

EqZ^' are identical in order of magnitude. Thus, the factor K of 

equation (5) for the formation of droplets appears smaller by l/n than 
for the formation of the crystal, where n denotes the molecule number 
of the nucleus. Although n is the order of magnitude of 100, this 
factor is not decisive in the problem involving the critical supersatura- 
tion. Moreover, it would be considerably overbalanced by the factor 


matter even less. 


omitted at ^crystal' factor A" and 

As long as no direct measurements of J are planned, but merely the order 
of magnitude of the critical 'supersaturation, the simple resilLt is : The 

factor K in Volmer's nucleation formula is simply equal to the manber 
of gas kinetic collisions, for the droplet as for the crystal. This 
statement applies, as seen in section 5 and section U, to linear and plane 
nuclei; naturally, involved here is solely the nvanber of collisions per 
second at the base. 


The decisive reason for the validity of the law of stages remains 

ctF 

then solely the fact that in the quantity rr— the surface F of the 

5kT 

nucleus corresponding to a certain super satiarat ion is greater on the cube 
than on the sphere. The difference in shape is the deciding factor, not 
the crystalline struct\xre. Its effect is computed on the assumption that 
the molecular volume v and the surface tension a for fluid and crystal 
are equivalent. 

If F = Cn2/3 is the surface corresponding to the molecule nimiber 

Pn 

n, then by equation (l), with x denoting the abbreviation of In — , 

Poo 


kTx 


^ = 2 oC^/2 
dn 5 pl/2 
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hence the surface corresponding to x 




9 (kTx)^ 

The nucleus volume V for the sphere (radius r) is 


V = nv = ~r 


5-, thus, F = 


for the cube (edge length a) : 


V 


= nv = a? y hence, F = 6v^/5n^/^ 


Hence, for the sphere 


c5 = 36rtv^ 


and for the cube 


C5 = 65y2 


The critical area corresponding to the same x is — = I. 9 I times 

greater for the cube than the droplet. To assure identical nucleation 
frequency, hence, equal values of F, it must 

. f f = /iTsi = 1.38 

^sphere \C sphere/ 

For the critical supersaturations themselves the condition would be 


JB- 

Poo 


cube 


‘^“^sphere 


( 58 ) 


As an illustration for applying this relationship, the supercooling 
at which crystalline and fluid nuclei occur with comparable frequency 
is analyzed. The saturation vapor pressure of the liquid phase is denoted 
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tiy Pi, that of the solid phase by P 2 - ^y equation ( 38 ) the condition 
for comparable nucleation frequency reads 

P2 \Pl/ 

or 

In p = In PI + 2.6(ln - In P 2 ) (39) 

(The factor 2.6 is equal to 1: (I .38 - 1). Figure 7 shows the vapor pres- 
sure curves In p^^ and In P 2 plotted against T. According to equa- 
tion ( 39 ) the curve for In p would then have about the shape of that 
indicated by the broken curve. Below this limiting curve, more crystal- 
line nuclei, above it, more fluid nuclei are to be expected. However, 
this theoretically interesting solution is meaningless in practice as 
long as the nucleation frequency lies below a limit amenable to observa- 
tion. For that reason it is necessary to determine, in the same manner 
as in section 2, the CTorve of that pressure at which a formation of fluid 
nuclei occurs at all in observable amounts (dotted curve). The intersec- 
tion point A of the two curves characterizes the temperature T^, at 
which an isothermal pressure rise would result in a simultaneous separa- 
tion of fluid and crystalline nuclei. Below T^ only solid, above T;^ 
only fluid nuclei would be observed. 

Naturally, it may also happen that no Intersection point appears. 

In that event, the law of stages holds unrestrictedly.® 


7 . THE GENERAL RESISTANCE ANALOGY 


As already stated several times in the foregoing, the equations (17)> 
( 23 ), and ( 32 ), applicable to the elementary process, can be so trans- 
formed by extension with a suitable factor that they could be interpreted 
as the Kirchhoff equations of a suitably chosen network of wires. It can 
be proved that this electrotechnical analogy is possible in complete 
generality for the condensation and dissolution process of any structure 
consisting of atoms. Again it is assumed that, besides the vapor phase 
of a substance, some fractions of another phase are present in a con- 
tainer. These fractions are hereinafter called crystals, without in any 
way infringing upon the general character. An uninterrupted input or 
transport of vapor and removal or addition of crystals of random specific 
size assures the steady distribution of the crystals of various sizes and 
shapes . 

®Such a case seems to exist in the theoretical case treated by 
Stranski and Totomanow (Z. f. phy. Chem. (A), I 63 , p. 399> 1933)- 
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Next, we consider any random specified type of crystals, say, of 
the shape represented in figure 8, for example. For a full description 
of such a crystal, a greater number of parameters are usually advantageous, 
a single one of which is, say, the number v of atoms in this crystal. 

By deposition of an atom at a well defined spot of this crystal, a crys- 
tal of type II with v + 1 atoms is produced. J is the excess per sec- 
ond of the growth processes which lead from I to II, through the evapora- 
tion processes which lead from II to I. Then, if nj and npp are 

the number of crystals of type I and II in the steady state, the equation 
for this specific transition process reads 


J = nja - npiq^+i 


( 1 * 0 ) 


where a and qy+1 are the repeatedly employed deposition and evapora- 
tion probabilities of the atom at that particular spot. With the abbre- 
viation Pv+i = — — the result is again 
^v+1 

i^II = i^lPv+l “ a Pv+1 

^V+l \|f 

Again qy+l = and a = F^ip^e ^ are introduced, hence, 

9y+i - t 

Pv+l = e M . 

The energy of separation cpy+l depends, as a rule, on all the param- 
eters of state of the states I and II, rather than on v alone. 

Now, visualize the crystal II built up successively from single 
atoms. To each one of the v + 1 single processes, there corresponds 
a specific Pj_. The individual still will be dependent upon the 

sequence in which the atoms of the crystal are Joined together. But the 
product P(ii) of all pp 
v+1 


v+1 

p(ii) = riPi = ® 

v+1 i=l 


^ 9 

^2^i-(v+l)t 

1=1 

kT 


^Factor p^^, which by itself corresponds to no growth process, is 
put equal to Ij hence = i]f. 
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is solely dependent upon the configiuration of crystal II, because the 

v+1 

total work of growth can no longer depend upon the manner in 

i=l 

which the growth took place. 

Dividing equation (4l) by Pjj gives 

v+1 


^I ~ ^II “ JRj (^2) 


where, for abbreviation 


and 




i -) O' 

'v(l) 


"II 

Py+ldl) 


(42a) 



(42b) 


Every possible form of crystal can now be characterized by a point 
in the space of the parameters, which define this form. To every possible 
transition, I — »II, we then correlate a wire connection between the 
points of state I and II, to which the resistance given by equation (4lb) 
is allocated. The equation (42) can then be regarded as the Kirchhoff 
equations of this wire netting and 0 as the corresponding potential of 
the nodes in this net. This interpretation is possible, because $ is 
merely dependent upon the state, but not the manner in which a crystal 
is built up. 

Obviously, this network of wires does not have to be multidimensional. 
Since the number of possible forms is finite, theoretically one parameter 
that counts the possible forms, may be sufficient. But for the represen- 
tation the use of two or three parameters, as in sections 4 and is more 
convenient so that the net becomes two or three dimensional. 

This network itself is rather complicated even in simple cases. In 
the networks treated in sections 4 and 5^ a large mmiber of wires were 
consistently ignored because of their high resistance and complete wire 
systems combined into one resistance. For the actual calculation, this 
general analogy is therefore of little help. However, the purely qijali- 
tative distribution of the resistances will be indicated. 



NACA 1374 


35 


The resistance of a wire depends solely upon its initial point and 
has the form 


V 


R 

V 

where the quantity vijf - 2 Z<Pi 

i=l 

formed to produce the system corresponding to the initial point of the 
wire by reversible process from the vapor. Of all the wire joints which 
lead from crystals with v atoms to those with v + 1 atoms, the wires 
proceeding from the crystals with the least work of growth are therefore 
the wires of the smallest resistance. Whether this minimum is always as 
sharp in more general cases as on the model used above, requires firrther 
study. During the advance along the road of minimum resistance from 
smaller to larger crystals, the work of growth must, at some time, reach 
a maximm value, because, while for very small systems it certainly 
increases with v, it must, at very great v become proportional with 
V negatively arbitrarily great, so far as the vapor is supersaturated 
at all with respect to the very large crystals. The crystal on which 
the resistance reaches its (absolute) maximum with regard to advancing 
with V and a minimum in comparison to the other wires with the same v, 
is called the Volmer nucleus. The resistance at this point is 

% 

^nucleus = | (^ 5 a) 

with Aj^ the work of nucleation. As is seen, the saddle-like character 
of the resistance distribution near the nucleus is completely independent 
of the model. The specific model representations merely yield information 
about the number of parallel wires of equal resistance in the saddle 
point, the distinct character of the saddle, and the extent to which any 
secondary maximums in the otherwise very jagged curve of the resistance 
become evident. The order of magnitude of the total resistance between 
two points with very small and very great v is, however, solely defined 

^nucleus . 


v\lr-^9i 

T i=l 

= - e kT (43) 

a 

in the exponent is the work to be per- 


Translated by J. Vanier 
National Advisory Committee 
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Figure 3.- A specific state of crystal growth. 
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Figure 4, ~ Network of partial currents for the growth of the linear chain. 
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Figure 5 .- 


Derivation of the relation PqP = 









state I 

a. 



State n 

d. 

Figure 8.- The most elementary process of crystal growth. Deposi- 
tion of one atom at the emphasized spot of state I leads to state n. 
Evaporation of raised atom on crystal II leads to crystal I. 
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